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Light-front quantum chromodynamics may lead to an accurate constituent ap- 
as . proximation for the low-energy properties of hadrons. This requires a cutoff that 

violates explicit gauge invariance and Lorentz covariance, leading to the calcula- 
OS ' tion of a renormalized QCD hamiltonian using a similarity renormalization group. 

Renormalization repairs broken symmetries, and in light- front field theory it moves 
dynamical effects that usually require large numbers of partons to few-body effec- 
tive interactions. This has been shown to work in QED through dominant contri- 
butions to the Lamb shift. In QCD logarithmic confinement arises at second order, 
and initial bound state calculations produce reasonable results. 
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1 Motivation and Outline 



■ Quantum chromodynamics (QCD) is the fundamental theory of the strong 

qq , interaction, but our understanding of QCD and our ability to use it to solve 

low-energy problems where the interactions are truly strong falls far short of 
our accomplishments in the study of quantum electrodynamics (QED). Both 
are gauge theories, but the fact that gluons carry color charge drastically 
complicates QCD. Couplings in QED are weak at low energies and photons 
are nearly free, so that the interactions their exchange produces are readily 
approximated. As a result low energy bound states can be accurately de- 
scribed using a small number of constituents whose interactions appear in the 
hamiltonian at second order, and the vacuum has no effect on bound states 
in light-front QED. 



The study of QCD and strong interaction phenomenology leads to little 
hope of deriving similar approximations. Interactions in QCD are strong at 
low energies and the color-charged gluons interact strongly while mediating in- 
teractions. The vacuum is supposed to have a complicated structure to which 
conferences are devoted. Its structure is assumed to be responsible for essen- 
tial aspects of the theory such as confinement and chiral symmetry breaking 
in all widely respected treatments of the theory, and it is hard to imagine 
an accurate description of individual hadrons that does not also include a 
complicated description of the vacuum in which they reside. Nonetheless, we 
advocate an approach in which an accurate .description of hadrons resembles 
the accurate description of atoms in QEDBEl 

We start with the heretical conjecture that a constituent picture of 
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hadrons can be derived from QCD. This conjecture guides our calculations, 
but the approach I describe is completely fixed by QCD and the constituent 
picture will fail if it is inadequate.^ 

If a constituent approximation is accurate, we can study the low-energy 
properties of hadrons (e.g., mesons) by solving a relativistic Schrodinger equa- 
tion: 

Ha I *a) = E | tf A ), (1) 

with, 

I *a) = <4 I qq) + < fg I qqg) +■•-, (2) 

where I use shorthand notation for the Fock space components of the state. 
The exact state vector includes an infinite number of terms. In a constituent 
approximation we truncate this series, adding terms to improve the approx- 
imation. We derive the hamiltonian from QCD, so we must allow for the 
possibility of constituent gluons. I have indicated that the hamiltonian and 
the state both depend on a cutoff, A, which is critical for the approximation. 

This approach has no chance of working without a renormalization scheme 
tailored to light-front hamiltonian field theory. Much of, our work has focused 
on the development of such a renormalization scheme.tli lj 

Consider the conditions under which it might be possible to truncate 
the above Fock space series without making an arbitrarily large error in the 
eigenvalue. I focus on the eigenvalue, because it is certainly not possible to 
approximate all observable properties of hadrons {e.g., wee parton structure 
functions) this way. For this approximation to be valid, all many-body states 
must approximately decouple from the dominant few-body components. 

We know that even in perturbation theory, high energy many-body states 
do not decouple from few-body states. In fact, the errors from simply discard- 
ing high energy states are infinite. In second-order perturbation theory, for 
example, high energy photons contribute an arbitrarily large shift to the mass 
of an electron. This second-order effect is illustrated in Fig. 1. The solution 
to this problem is well-known, renormalization. Renormalization moves the 
effects of high energy components in the state to effective interactions in the 
hamiltonian. 

It is difficult to see how a constituent approximation can emerge in a 
hamiltonian calculation using any regularization scheme without a cutoff that 

a The description of light hadrons requires refinements to the simple approximations I de- 
scribe, because of chiral symmetry breaking. 
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Figure 1. The second-order shift in the self-energy of a bare electron due to mixing with 
electron-photon states. 



either removes high energy degrees of freedom or removes direct couplings 
between low and high energy degrees of freedom. 

In the best case scenario we expect the cutoff to act like a resolution. If 
the cutoff is increased to an arbitrarily large value, the resolution increases 
and instead of seeing a few constituents we resolve the substructure of the 
constituents and the few-body approximation breaks down. As the cutoff is 
lowered, this substructure is removed from the state vectors, and the renor- 
malization procedure replaces it with effective interactions in the hamiltonian. 
Any "cutoff" that does not remove this substructure from the states is of no 
use to us. 

This._point is well-illustrated by the QED calculations discussed 
below.otffiJ There is a window into which the cutoff must be lowered for the 
constituent approximation to work. If the cutoff is raised atomic states start 
to include photons, and as the cutoff is raised further they start to include 
additional photons and electron-positron pairs. After the cutoff is lowered 
to a value that can be self-consistently determined a-posteriori, photons and 
pairs are removed from the states and replaced by the Coulomb interaction 
and relativistic corrections in the hamiltonian. The cutoff cannot be lowered 
too far using a perturbative renormalization group, hence the window. 

Thus, if we remove high energy degrees of freedom, or coupling to high 
energy degrees of freedom, we should encounter self-energy corrections leading 
to effective one-body operators, vertex corrections leading to effective vertices, 
and exchange effects leading to explicit many-body interactions not found in 
the canonical hamiltonian. We naively expect these operators to be local 
when acting on low energy states, because simple uncertainty principle argu- 
ments indicate that high energy virtual particles cannot propagate very far. 
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Figure 2. Mixing of low-energy few-body states with high-energy many-body states alters 
the dispersion relation for single particles, alters emission and absorption vertices, and 
produces approximately local few-body interactions. 



Unfortunately these arguments break down in light-front coordinates, and at 
best we can maintain transverse locality.0 

In this article I give a brief overview of the new renormalization techniques 
we employ, and in a separate article in these Proceedings liil Sergio Szpigel 
and I illustrate them using the two-dimensional <5-function potential. 

Low energy many-body states do not typically decouple from low energy 
few-body states. The worst of these low energy many-bodi. states is the 
vacuum. This is what drives us to use light-front coordinates.E3 Fig. 3 shows 
a pair of particles being produced out of the vacuum in equal-time coordinates 
t and z. The transverse components x and y are not shown, because they are 
the same in equal-time and light-front coordinates. The figure also shows the 
light-front time axis, 

x + = t + z , (3) 

and the light-front longitudinal spatial axis, 

x~ = t - z . (4) 

In equal-time coordinates it is kinematically possible for virtual pairs to 
be produced from the vacuum, as long as their momenta sum to zero so that 
three-momentum is conserved. Because of this, the state vector for a proton 
includes an arbitrarily large number of particles that are disconnected from 
the proton. The only constraint imposed by relativity is that particle velocities 
be less than or equal to that of light. 
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Figure 3. Light-front coordinates. Light-front 'time' is x + = t + z, and light-front longitu- 
dinal position is a: - = t — z. 



In light-front coordinates, however, we see that all allowed trajectories lie 
in the first quadrant. In other words, light-front longitudinal momentum, p + 
(conjugate to x~ since a ■ b — ±(a + fr~ + a~b + ) — aj^ • bj_), is always positive, 

P + > . (5) 

We exclude particle modes with p + = 0, forcing the vacuum to be trivial 
because it is the only state with p + = 0. Moreover, the light-front energy of 
a free particle of mass m is 

- pj + m 2 

This implies that all free particles with zero longitudinal momentum have 
infinite energy, unless their mass and transverse momentum are zero. 

Is the vacuum really trivial? What about confinement? What about 
chiral symmetry breaking? What about instantons? What about the job 
security of theorists who study the vacuum? 

I simply discard all p + — degrees of freedom and replace their effects 
using a renormalization procedure that does not require any explicit reference 
to zero modes. Thus the vacuum in our formalism is trivial. We are forced 
to work in the "hidden symmetry phase" of the theory, and to introduce 
effective interactiQna.±hat reproduce all effects associated with the vacuum in 
other formalisms. EltJO The simplest example of this approach is provided by 
a scalar field theory with spontaneous symmetry breaking. It is possible to 
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shift the scalar field and deal explicitly with a theory containing symmetry 
breaking interactions. In the simplest case <f> 3 is the only relevant or marginal 
symmetry breaking interaction, and one can simply tune this coupling to the 
value corresponding to spontaneous rather than explicit symmetry breaking. 

The use of a symmetry-breaking cutoff and the removal of zero-modes 
leads to a large number of operators not found in the canonical QCD hamil- 
tonian. This complicates the renormalization procedure, but it may lead to 
tremendous simplifications in the final nonperturbative problem. For exam- 
ple, few-body operators must produce confinement manifestly! 

Confinement cannot require particle creation and annihilation, flux tubes, 
etc. This is easily seen using a variational argument. Consider a color neutral 
quark-antiquark pair that are separated by a distance R, which is slowly 
increased to infinity. Moreover, to see the simplest form of confinement assume 
that there are no light quarks, so that the energy should increase indefinitely 
as they are separated if the theory possesses confinement. At each separation 
the gluon components of the state adjust themselves to minimize the energy. 
But this means that the expectation value of the hamiltonian for a state 
with no gluons must exceed the energy of the state with gluons, and therefore 
must diverge even more rapidly than the energy of the true ground state. This 
means that there must be a two-body confining interaction in the hamiltonian. 
If the renormalization procedure is unable to produce such confining two-body 
interactions, it is invalid. 

1.1 Simple Strategy 

I want to outline a conceptually simple strategy for bound state calculations S 
The first step is to use. a perturbative similarity renormalization group Em and 
coupling coherence ErtH to find the renormalized hamiltonian as an expansion 
in powers of the canonical coupling: 



We compute this series to a finite order and to date have not required any 
ad hoc assumptions to uniquely fix the hamiltonian. No operators are added 
to the hamiltonian by hand, so it is completely determined by the underlying 
theory to this, order. This step is illustrated in a separate article in these 



The second step is to employ bound state perturbation theory to solve 
the eigenvalue problem. The complete hamiltonian contains every interaction 
(although each is cut off) contained in the canonical hamiltonian, and many 



H A = h + g A h A + g\h$ + ■■■. 



(7) 
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more. We separate the hamiltonian, 

H A = W A + V A , (8) 



treating 7i A nonperturbatively and computing the effects of V A in bound 



state perturbation theory. We must choose Hq and A so that Hq is tractable 
and to minimize corrections from higher orders of V A within a constituent 
approximation. 

If a constituent approximation is valid after A is lowered to a critical value 
that must be determined, we may be able to move all particle creation and 
annihilation to V A . Hq includes many-body interactions that do not change 
particle number, and these interactions should be primarily responsible for 
the constituent bound state structure. 

There are several obvious flaws in this strategy. Chiral symmetry- 
breaking operators, which must be included in the hamiltonian since we work 
entirely in the hidden symmetry phase of the theory, do not appear at any fi- 
nite order in the coupling. There is only one relevant chiral symmetry breaking 
operator, and it appears in the canonical hamiltonian when quarks are mas- 
sive (spin-flip gluon emission by quarks) although it can acquire non-canonical 
dependence on longitudinal momenta since there is no longitudinal locality. 
This operator must simply be added if quarks are massless and tuuecLia fit 
spectra or fixed by a non-perturbative renormalization procedure. EloLZI In 
addition, there are perturbative errora-in the strengths of all operators. We 
know from simple scaling arguments £3 that when A is in the perturbative 
scaling regime: 

• small errors in relevant operators exponentiate in the output, 

• small errors in marginal operators produce comparable errors in output, 

• small errors in irrelevant operators tend to decrease exponentially in the 
output. 

This means that even if a relevant operator appears (e.g., a constituent quark 
or gluon mass operator), we may need to tune its strength rather than use 
its perturbative value to obtain reasonable results. We have not had to do 
this, but pjze have recently studied some of the effects of tuning a gluon mass 
operatorEl 

To date this strategy has produced well-known results in QED BEEj 
through the Lamb shift, and reasonable results for heavy quark bound states 
in QCD.ElBtj All of these calculations rely on a nonrelativistic reduction 
of the effective hamiltonian, which leads to drastic simplifications. Glueball 
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calculations are being completed 123 and will be the first fully relativistic QCD 
calculations in this approach. 

The best place to begin the study of light-front field theory and to see 
how a constituent approximation can arise in light-fron±_gauge theories is the 
Schwinger model, massless QED in 1 + 1 dimensions.^ This model can be 
solved analytically] 2 -! and its physical content is remarkably simple. There 
is only one physical particle, a massive neutral scalar particle with no self- 
interactions. The Fock space content of the physical states depends crucially 
on the coordinate system and gauge, and it-is. only in light-front coordinates 
that a simple constituent picture emerges.E3o In light-front field theory the 
physical particle is a bound state of a single electron-positron pair, with a 
wave function that is constant in the longitudinal fraction carried by either 
particle. 

The Schwinger model does not display the renormalization problems that 
must be solved in QED3+1 and QCD3+1. Before turning to these theories, I 
discuss the renormalization machinery that has been developed for light-front 
hamiltonians. 

2 Light-Front Renormalization Group 

In 3 + 1 dimensions we must introduce a cutoff, A, and we never perform 
explicit bound state calculations with A anywhere near its continuum limit. In 
fact, we want to let A become as small as possible. In my opinion, any strategy 
for solving light-front QCD that requires the cutoff to explicitly approach 
infinity in the nonperturbative part of the calculation is useless. Therefore, 
we must set up and solve 

P 2 4- M 2 

PK I *a(^)> = p + I *a(P)> • (9) 

Physical results, such as the mass, M, can not depend on the arbitrary 
cutoff, A. This means that P^ and | \&a) must depend on the cutoff in such 
a way that (\&a I P\ I ^a) does not. Wilson b a i sprl the derivation of his 
renormalization group on this observation,L3EZlt3lla and we modify Wilson's 
renormalization group to compute . 

It is difficult to even talk about how the hamiltonian depends on the 
cutoff without having a means of changing the cutoff. If we can change the 
cutoff, we can explicitly watch the hamiltonian's cutoff dependence change 
and fix its cutoff dependence by insisting that this change satisfy certain 
requirements {e.g., that the limit in which the cutoff is taken to infinity exists). 
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Figure 4. Two ways to run a cutoff on free energy. In (a) a cutoff on the magnitude of 
the energy is lowered from the solid to the dashed lines, with problems resulting from the 
removed shaded region. In (b) a cutoff on how far off diagonal matrix elements appear is 
lowered from the dashed to the solid lines. 



We introduce an operator that changes the cutoff, 



where I assume that Ai < Ao. To simplify the notation, I let H(Ai) — Hi. To 
renormalize the hamiltonian we study the properties of the transformation. 

Fig. 4 displays two generic cutoffs that might be used. Traditionally the- 
orists have used cutoffs that remove high energy states, as shojsai in Figure 
4a. This is the type of cutoff Wilson employed in his initial work c3 and I have 
studied its use in light-front field theory.D When a cutoff on energies is reduced, 
all effects of couplings eliminated must be moved to effective operators. When 
these effective operators are computed perturbatively they involve products of 
matrix elements divided by energy denominators. Expressions closely resem- 
ble those encountered in standard perturbation theory, with the second-order 
operator involving terms of the form 



This new effective interaction replaces missing couplings, so the states 4>i and 
4>j are retained and the state 4>k is one of the states removed. The problem 
comes from the shaded, lower right-hand corner of the matrix, where the 
energy denominator vanishes for states at the corner of the remaining matrix. 
In this corner we should use nearly degenerate perturbation theory rather than 
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H(A0 = T[H(A )} , 



(10) 




(ii) 



perturbation theory, but to do this requires solving high energy many-body 
problems nonperturbatively before solving the low energy few-body problems. 

An alternative cutoff, which does not actually remove any states and that 
can be run by a similarity transformation^ is shown in Fig. 4b. This cutoff 
removes couplings between states whose free energy differs by more than the 
cutoff. The advantage of this cutoff is that the effective operators resulting 
from it contain energy denominators that are never smaller than the cutoff, 
so that a perturbative approximation for the effective hamiltonian may work 
well. 

Before discussing the similarity transformation, which runs a cutoff on 
off-diagonal matrix elements, I need to introduce the classification of opera- 
tors as relevant, marginal, and irrelevant. This classification is not necessary 
to use the similarity renormalization group, but it is necessary for the defini- 
tion of coupling coherence, and it is essential for a full understanding of the 
renormalization group. It starts with the definition of a fixed point, which is 
a hamiltonian that is left invariant by the transformation, 

H*=T[H*\. (12) 

Obviously a fixed point is a renormalized hamiltonian, because the cutoff 
fixed point hamiltonian is identical to the fixed point hamiltonian with infinite 
cutoff. 

Consider the immediate 'neighborhood' of the fixed point, and assume 
that the trajectory remains in this neighborhood. This assumption must be 
justified a posteriori, but if it is true we should write 

ff, + (13) 

and consider the trajectory of small deviations 5 Hi. 

As long as SHi is 'sufficiently small,' we can use a perturbative expansion 
in powers of 8 Hi, which leads us to consider 

6Hi +1 = L ■ SHi + N[6Hi] . (14) 

Here L is the linear approximation of the full transformation in the neigh- 
borhood of the fixed point, and N[SHi] contains all contributions to SHi + i of 
0{5H?) and higher. 

The object of the renormalization group calculation is to compute trajec- 
tories and this requires a representation for 5 Hi. The problem of computing 
trajectories is one of the most common in physics, and a convenient basis 

'in deference to the original work I call this a similarity transformation even though in all 
cases of interest to us it is a unitary transformation. 
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for the representation of SHi is provided by the eigenoperators of L, since L 
dominates the transformation near the fixed point. These eigenoperators and 
their eigenvalues are found by solving 

L ■ O m — \ m O m ■ (15) 

For light-front field theory this linear transformation is a scaling of the trans- 
verse coordinate, the eigenoperators are products of field operators and trans- 
verse derivatives, and the eigenvalues are determined by the transverse di- 
mension of the operator. All operators can include both powers and inverse 
powers of longitudinal derivatives because there is no longitudinal locality.!!] 
Using the eigenoperators of L as a basis we can represent SHi , 

SHi = ^ VrniOm + ^ 9r ni O m + ^ w m t O m . (16) 

roei? mEM mel 

Here the operators O m with m € R are relevant (i.e., X m > 1), the opera- 
tors O m with m S M are marginal (i.e., X m — 1), and the operators with 
to € I are irrelevant (i.e., X m < 1). The motivation behind this nomencla- 
ture is made clear by considering repeated application of L, which causes the 
relevant operators to grow exponentially, the marginal operators to remain 
unchanged in strength, and the irrelevant operators to decrease in magnitude 
exponentially 



2.1 Similarity transformation 

Stan Glazek and Ken Wilson studied the problem oLamall energy denomi- 
nators that typically appear in effective interactions,^ and realized that a 
similarity transformation that runs a different forn^of cutoff (as discussed 
above) avoids this problemB Independently, Wegnera developed a similarity 
transformation that is easier to use than those first studied by Glazek and 
Wilson. Recently Walhout has developed a transformation that may be well- 
suited to analytical and higher-order calculations.^ 

Many details and a simple example of the use of Wegner's similarity trans- 
formation can be found elsewhere in these Proceedings .til Here I summarize 
the details needed to complete QED and QCD calculations through second 
order. Consider a hamiltonian, H s = h + V s , where h is a free hamiltonian 
that may contain mass terms. The dependence of H on the cutoff, where s is 
the inverse of the cutoff squared, is given by 

^ = [H s ,[H s ,h}. (17) 
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This equation is most easily analyzed in terms of matrix elements between 
eigenstates of h, 

h\<f>i) = 6i|0i) . (18) 

A reduced interaction, V a , can be denned, 

V sij = e~ sA *> V sij , (19) 

where Ay = Ci—ej. A gaussian cutoff factor that clearly forces the interactions 
towards the diagonal is isolated, and the equation for V s starts at second 
order, 

+ A* <») 

k 

where I use A 2 j - A 2 ik - A 2 jk = -2A ik A jk . 

This is a first-order differential equation and it is exact. Its solution 
is completely determined when a complete set of boundary conditions are 
specified. In its simplest form the boundary conditions would be given by 
specifying H So using a single value of sq for all matrix elements, with sq = 
corresponding to an infinite cutoff. As discussed in the simple 5-function 
example,EJ divergences prevent us from using sq = to specify boundary 
conditions for relevant and marginal operators, but we are free to choose 
different values of sq for different operators and this freedom is crucial. The 
boundary conditions are given by coupling coherence, which 1 discuss in the 
next section. 

For QED and QCD we can approximate V s using a perturbative expan- 
sion in powers of a single running coupling constant, 

_ _(1) „_(2) 

V s = g s V +g 2 s V s +-... (21) 

In both QED and QCD the coupling does not begin to run until third order, 

<i) 

and the first order term, V , is given by the canonical hamiltonian and has 
no dependence on s other than that of the running coupling. To second order 
we have 

(2) 

= E (A ik + A jk ) v2 <• e- 2 ^^ . (22) 

k 



Integrating from so to s, we obtain 



-< 2 > 1 v-^— W —CD / 1 1 , 



x [ e ~2s A lk A jk _ e - 2s A lfe A jfc j _ (23) 
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This simple equation is sufficient to compute the renormalized QED and QCD 
hamiltonians through second order the canonical couplings. 

In order to use the concept of fixed points in the similarity renormalization 
group, we need to modify the transformation slightly so that it is possible to 
get a fixed point with interactions. The free hamiltonian h is obviously a fixed 
point, but Eq. ( |20| ) clearly shows that the linearized transformation is given 
by Eq. (|l9|). If we want to have fixed points for which V s is not zero, we 
need to modify the transformation so that the gaussian cutoff factor in Eq. 
( |l9| ) does not change as s changes. This is most easily accomplished in the 
traditional manner by rescaling variables so that the eigenvalues in Eq. ( |l9| ) 
absorb the change in s. For example, in massless light-front field theory, 

= % ( 24 ) 

Pi 

We can absorb a change in s by rescaling all transverse momenta, and this 
leads to an operator classification based on the transverse engineering dimen- 
sion of operators. There are other interesting possibilities,!] but the result is 
that transverse local operators (i.e., no inverse powers of transverse momenta) 
that vanish when all transverse momenta are taken to zero are irrelevant. The 
electron-photon and quark-gluon couplings are marginal, excentibr a spin-flip 
piece that is relevant and breaks light-front chiral symmetry.EHl3 

The interesting new feature of light-front field theory that is not encoun- 
tered in equal-time or euclidean field theory is that there is no longitudi- 
nal locality, so that the longitudinal dimension of operators does not affect 
their classification. A relevant or marginal operator can contain a function 
of longitudinal momentum fractions, and, these functions inevitably appear 
in renormalized light-front hamiltonians .0 This means that there are effec- 
tively an infinite number of relevant and marginal operators in renormalized 
light-front hamiltonians, and such a situation is usually regarded as a disaster 
for renormalization group treatments because it indicates that there may be 
an infinite number of free physical parameters-Ken Wilson and I developed 
coupling coherence to deal with this problem.tiHH 

2.2 Coupling coherence 

The basic mathematical idea behindpCpupling coherence was first formulated 
by Oehme, Sibold, and Zimmerman.L3 They were interested in field theories 
where many couplings appear, such as the standard model, and wanted to 
find some means of reducing the number of couplings. 

The puzzle that led us to the same results is how to reconcile our knowl- 
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edge from covariant formulations of QCD that only one running coupling 
constant characterizes the renormalized theory with the appearance of new 
counterterms and functions required by the light-front formulation. What 
happens in perturbation theory when there are effectively an infinite number 
of relevant and marginal operators? In particular, does the solution of the 
perturbative renormalization group equations require an infinite number of 
independent counterterms (i.e., independent functions of the cutoff)? Cou- 
pling coherence provides the conditions under which a finite number of running 
variables determines the renormalization group trajectory of the renormalized 
hamiltonian. To leading nontrivial orders these conditions are satisfied by the 
counterterms introduced to restore Lorentz covariance in scalar field theory 
and gauge invariance in light-front gauge theories. In fact, the conditions 
can apparently be used to determine all counterterms in the hamiltonian, in- 
cluding relevant and marginal operators that contain functions of longitudinal 
momentum fractions; and with no direct reference to Lorentz covariance, this 
symmetry seems to be restored to observables.LJ 

A coupling-coherent hamiltonian is analogous to a fixed point hamilto- 
nian, but instead of reproducing itself exactly it reproduces itself in form with 
a limited number of independent running couplings. If g\ is the only inde- 
pendent coupling in a theory, in a coupling-coherent hamiltonian all other 
couplings are invariant functions of g\, fi(g\)- The extra couplings /j(<7a) 
depend on the cutoff only through their dependence on the running coupling 
g\, and in general we demand fi(0) = 0. This boundary condition on the 
dependent couplings is motivated in our calculations by the fact that it is the 
combination of the cutoff and the interactions that force us to add the coun- 
terterms we seek, so the counterterms should vanish when the interactions are 
turned off. 

Let me illustrate the idea with a simple example with a finite number 
of relevant and marginal operators ab initio, and use coupling coherence to 
discover when only one or two of these may independently run with the cutoff. 
Such conditions are met when an underlying symmetry exists, although this 
is not necessary. 

Consider a theory in which two scalar fields interact, 

V&) = ^<t>t + ±4 + ■ (25) 

Under what conditions are there fewer than three independent running cou- 
pling constants? We can use a simple cutoff on Euclidean momenta, q 2 < A 2 . 
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Letting t = ln(A/Ao), the Gell-Mann-Low equations are 
d\ 1 

^=3CA? + -CA^+0(2 loop), 
<9A 1 

^=3CA^ + -CA^+0(2 loop), 

-gf = + CAiA 3 + CA 2 A 3 + 0(2 loop) ; 

where £ = 7i/(167r 2 ). It is not important at this point to understand how 
these equations are derived. 

Assume that there is only one independent variable, A = Ai, so that A2 
and A3 are functions of A. In this case Eqs. ( p7|) and (pq ) become, 

d\ 2 



(26) 
(27) 
(28) 



3A 2 + ^A 2 



1 

12' 



3A 2 . + —A 



12 



^3 j 



(29) 



f3A 2 + -A 2 ) = I 
V 12 3 ^ £>A 3 



— r As + AA 3 + A2A3 . 



(30) 



The only non-trivial solutions are A2 = A, and either A3 = 2A or A3 = 6A. If 
A3 = 2A, 



A 
4! 



and we find the (9(2) symmetric theory. If A3 = 6A, 



v(4) = 



A 



2-4! 



(31) 



(32) 



and we find two decoupled scalar fields. Therefore, A2 and A3 do not run 
independently with the cutoff if there is a symmetry that relates their strength 
to Ai. 

The condition that a limited number of variables run with the cutoff does 
not only reveal symmetries broken by the regulator, it may also be used to un- 
cover symmetries that are broken by the vacuum. For example, it is straight- 
forward to show that in a scalar theory with a (j) 3 coupling, this coupling can 
be fixed as a function of the cf> 2 and </> 4 couplings _pnly if the symmetry is 
spontaneously broken rather than explicitly brokcn.tj 

For the QED and QCD calculations, I need to compute the hamiltonian 
to second order, while the canonical coupling runs at third order. In this 
case we can use Eq.(E3J), with sq — ► 00 for relevant operators and sq — > for 
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irrelevant operators, and with the bare coupling e replaced by the running 
coupling e s or e^y. The more interesting case of marginal operators can be 
avoided at first .LrtJ 



3 Light-Front QED and QCD 

Various forms of the canonical light-front QED and QCD hamiltonians can be 
found in several articles EilH EaLLowing Brodsky and Lepage, I have displayed 
these hamiltonians elsewhere .0H 

In light-cone gauge and using light-front coordinates, it is possible to ex- 
plicitly eliminate all unphysical degrees of freedom and write the hamiltonian 
in terms of two-component fermions and transverse gluons. Any ambiguities 
in the procedure that come from the zero-mode problem or normal-ordering 
are resolved by coupling coherence, so the renormalized hamiltonian is appar- 
ently uniquely determined order-by-order in the running coupling. 

3.1 Light-front QED 

In this section I follow the strategy outlined in the first section to compute 
the positronium spectrum. I outline the calculation through the leading order 
Bohr results B and indicate how higher order calculations proceedBtd 

The first step is to compute a renormalized cutoff hamiltonian as a power 
series in the running coupling eA, 

Hjtf =h + etJix + e\h 2 H \-e%h N . (33) 

Having obtained the hamiltonian to some order in e\, the next step is to split 
it into two parts, 

H A = H$ + V A . (34) 

Hq must be accurately solved non-perturbatively, producing a zeroth order 
approximation for the eigenvalues and eigenstates. The greatest ambiguities 
in the calculation appear in the choice of 7i A , which requires one of science's 
most powerful computational tools, trial and error. 

In QED and QCD I conjecture that for sufficiently small A all interactions 
in Hq preserve particle number, with all interactions that involve particle 
creation and annihilation in V A . Corrections from V A are then computed in 
bound state perturbation theory. 

Since 7i A is assumed to include interactions that preserve particle number, 
the zeroth order positronium ground state is a pure electron-positron state. 
We only need one- and two-body interactions; i.e., the electron self-energy 
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and the electron-positron interaction. The hamiltonian is computed to second 
order using Eq. (|2|). We must specify sq in Eq. (|23|), which corresponds to 
the inverse cutoff squared at which boundary conditions are placed on the 
hamiltonian. Coupling coherence leads to the prescription that sq — ► for 
irrelevant operators and sq — ► oo for relevant operators. 

Bare electron mixing with electron-photon states leads to a self-energy 
(see Fig. 1): 



a e 2 A2 / y 2 A 2 \ 3 1 ym 2 A 2 

Zih(p) = < 2yA 2 In - - y A 2 



8ir 2 p+ 1 I (yA 2 + m 2 )e/ 2 2 yA 2 + m 2 

+ 3 TO 2 ln(-^)} + (e / y ); (35) 

where y is the fraction of longitudinal momentum carried by the electron, 
y = p + 1 P + . To simplify the discussion I have replaced the gaussian cutoff 
factors that appear in all integrals with step functions, and completed the 
integrals analytically using l/*/s = A 2 /P + . It isjpossible to produce such step 
function cutoffs with a similarity transformation,Elbut this leads to pathologies 
at higher order. More importantly, I have been forced to introduce a second 
cutoff, 

xp+ > eP + , (36) 

because there is a logarithmic divergence in the loop longitudinal momentum 
integration even with the gaussian cutoff in place. This second cutoff must 
be taken to zero and no new counterterms can be added to the hamiltonian, 
so all divergences must cancel before it is taken to zero. 

We have no choice about whether this divergent operator is in the hamil- 
tonian if we use coupling coherence. We can only choose between putting it in 
H.q or in V A . I make different choices in QED and QCD, and the arguments 
are based on physics. 

The divergent electron 'mass' is a complete lie. We encounter a term 
proportional to e^A 2 ln(l/e)/P + when the scale is A; however, we can reduce 
this scale as far as we please in perturbation theory. Photons are massless, so 
the electron continues to dress itself with small-x photons to arbitrarily small 
A. Since I believe that this divergent self-energy is exactly canceled by mixing 
with small-x photons, and that this mixing is perturbative in QED, I simply 
put it in V A . 

There are two time-ordered diagrams involving photon exchange between 
an electron with initial momentum p± and final momentum p2 , and a positron 
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Figure 5. Effective two-body interaction between low energy constituents resulting from: (i) 
the canonical instantaneous exchange interaction, and (ii) the elimination of direct coupling 
between low energy two-body states and high energy states containing an additional gauge 
particle. 



with initial momentum ki and final momentum ki- These are shown in Fig. 
5, along with the instantaneous-exchange diagram. I refer the reader to longer 
articles where details are given 013 and concentrate here on the essential results. 

Photon exchange above the cutoff produces an effective interaction that 
cancels the instantaneous photon exchange operator in the canonical hamilto- 
nian, replacing it with a Coulomb interaction and magnetic interactions that 
are partially responsible for fine structure. Instantaneous photon exchange 
below the cutoff remains, and as is discussed in the section on QCD it pro- 
duces a long-range confining interaction. In QED this long range interaction 
is exactly cancelled by further exchange of massless low energy photons in 
bound state perturbation theory. In QCD the long range interaction also acts 
on gluons, blocking such a cancellation and producing the essential difference 
between QED and QCD in this approach. 

This means that we can concentrate on photon exchange above the cutoff 
to leading order, which still leaves us with a complicated operator in light- 
front coordinates. In order to present an analytic analysis I make assumptions 
that are justified a posteriori. First I assume that the electron and positron 
momenta can be arbitrarily large, but that in low-lying states their relative 
momenta satisfy 

|p_L-k x | ~am, (37) 
\ P + -k+\~a(p+ + k+) . (38) 
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This allows us to use power counting to evaluate the perturbative strength of 
operators for small coupling. 

Given these order of magnitude estimates for momenta, we can drastically 
simplify all of the operators in the hamiltonian. At this point we can complete 
the zeroth order analysis of positronium using the state, 

- e / tffi tjz (P - v - » 

0(p,a;k,\)b\p,<j)J(k,\)\O) , (39) 



where <p(p, a; k, A) is the wave function for the relative motion of the electron 
and positron, with the center-of-mass momentum being P. We need to choose 
the longitudinal momentum appearing in the cutoff, and I use the natural scale 
P+. 

If we want to find a cutoff for which the ground state is dominated by 
the electron-positron component of the wave function, we need the cutoff to 
remove the important part of the electron-positron-photon phase space. Since 
the exchanged photon energy is typically 0(am 2 ), we need: 

A 2 < am 2 . (40) 

On the other hand, wc cannot allow the cutoff to remove the region of the 
electron-positron phase space from which the wave function receives most of 
its strength. This requires 

A 2 > a 2 m 2 . (41) 

For cutoffs that satisfy am 2 > A 2 > a 2 m 2 , the bound state equation can 
be simplified to: 

k 2 f r1 3 ki 1 

- E ^ = ^ - a J R(i^)^ (k2) ' (42) 

where a simple change of variables has replaced longitudinal momentum frac- 
tions with a z-component of momentum, making the system's nonrelativistic 
dynamics manifest. The cutoffs drop out to leading order, leaving us with the 
familiar nonrelativistic Schrodinger equation for positronium in momentum 
space. The solution is 

<Hk) = N , (43) 
(k 2 +mE) 

E = -a 2 m . (44) 
4 
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AT is a normalization constant. This is the Bohr energy for the ground state 
of positronium, and it is obvious that the entire nonrelativistic spectrum is 
reproduced to leading order. 

Beyond this leading order result the calculations become much more inter- 
esting, and in any hamiltonian formulation they rapidly become complicated. 
The leading correction to the binding energy is 0(a 4 ), and producing these 
corrections is a much more serious test of the renormalization procedure. We 
have shown that the fine structure of positronium is correctly reproduced when 
the first- and .second-order corrections from bound state perturbation the- 
ory are added.y This is a formidable calculation, because the exact Coulomb 
bound and scattering states appear in second-order bound state perturbation 
theory. 

A complete calculation of the Lamb shift in hydrogen would also require 
a fourth-order similarity calculation of the hamiltonian; however, the dorxtj 
inant contribution to the Lamb shift that was first computed by Bethe 
can be computed using a hamiltonian determined to C(a).E3 In our calcula- 
tion a Bloch transformation was used rather than a similarity transformation 
because the Bloch transformation is simpler and small energy denominator 
problems can be avoided in analytical QED calculations. 

The primary obstacle to using our light-front strategy for precision QED 
calculations is algebraic complexity. We have successfully used QED as a 
testing ground for this strategy, but these calculations can be done much 
more conveniently using other methods. The theory for which we believe our 
methods are best suited is QCD. 

3.2 Light-front QCD 

We only require the QCD hamiltonian determined to 0(a) to discuss a simple 
confinement mechanism that appears naturally in light-front QCD and to 
complete reasonable zeroth order calculations for heavy quark bound states. 
To this order the QCD hamiltonian in the quark- antiquark sector is almost 
identical to the QED hamiltonian in the electron-positron sector. Of course 
the QCD hamiltonian differs significantly from the QED hamiltonian in other 
sectors, and this is essential for justifying my choice of TLq for non-perturbative 
calculations. 

The basic strategy for doing a sequence of (hopefully) increasingly accu- 
rate QCD bound state calculations is almost identical to the strategy for doing 
QED calculations. Find an expansion for H A in powers of the QCD coupling 
constant to a finite order. Divide the hamiltonian into a non-perturbative 
part, H.Q , and a perturbative part, V A . The division is based on the physical 
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argument that adding a parton in an intermediate state should require more 
energy than indicated by the free hamiltonian, and that as a result these states 
'freeze out' as the cutoff approaches Aqcd- When this happens the evolution 
of the hamiltonian as the cutoff is lowered further changes qualitatively, and 
operators that were consistently canceled over an infinite number of scales also 
freeze, so that their effects in the few parton sectors can be studied directly. 
A one-body operator and a two-body operator arise in this fashion, and serve 
to confine both quarks and gluons. 

The simple confinement mechanism I outline is certainly not the final 
story, but it may be the seed for the full confinement mechanism. One of the 
most serious problems we face when looking for non-perturbative effects such 
as confinement is that the search itself depends on the effect. A candidate 
mechanism must be found and then shown to produce itself self-consistently 
as the cutoff is lowered towards Aqcd- 

Once we find a candidate confinement mechanism, it is possible to study 
heavy quark bound states with little modification of the QED strategy. Of 
course the results in QCD differ from those in QED because of the new choice 
of Hq , and in higher orders because of the gluon interactions. 

When we compute the QCD hamiltonian to 0(a), several significant new 
features appear. First are the familiar gluon interactions. In addition to the 
many gluon interactions found in the canonical hamiltonian, there are modi- 
fications to the instantaneous gluon exchange interactions, just as there were 
modifications to the electron-positron interaction. For example, a Coulomb 
interaction automatically arises at short distances. In addition the gluon self- 
energy differs drastically from the photon self-energy. 

The photon develops a self-energy because it mixes with electron-positron 
pairs, and this self energy is 0{aK 2 / P + ). When the cutoff is lowered below 
4m 2 , this mass term dies exponentially because it is no longer possible to 
produce electron-positron pairs. For all cutoffs the small bare photon self- 
energy is exactly canceled by mixing with pairs below the cutoff. I do not go 
through the calculation, but because the gluon also mixes with gluon pairs 
in QCD, the gluon self-energy acquires an infrared divergence, just as the 
electron did in QED. In QCD both the quark and gluon self-energies are 
proportional to aA 2 ln(l/e)/P + , where e is the secondary cutoff on parton 
longitudinal momenta introduced in the last section. This means that even 
when the primary cutoff A 2 is finite, the energy of a single quark or a single 
gluon is infinite, because we are supposed to let e — > 0. 

In QED I argued that the bare electron self-energy is a complete lie, 
because the bare electron mixes with photons carrying arbitrarily small lon- 
gitudinal momenta to cancel this bare self-energy and produce a finite mass 
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physical electron. However, in QCD there is no reason to believe that this 
perturbative mixing continues to arbitrarily small cutoffs. There are no mass- 
less gluons in the world. In this case, the free QCD hamiltonian is a complete 
lie and cannot be trusted at low energies. 

On the other hand, coupling coherence gives us no choice about the quark 
and gluon self-energies as computed in perturbation theory. The question is 
not whether large self-energies appear in the hamiltonian. The question is 
whether these self-energies are canceled by mixing with low energy multi-gluon 
states. As the cutoff approaches Aqcd, I speculate that these cancellations 
cease to occur because perturbation theory breaks down and a mass gap 
between states with and without extra gluons appears. 

But if the quark and gluon self-energies diverge, and the divergences can- 
not be canceled by mixing between sectors with an increasingly large number 
of partons, how is it possible to obtain finite mass hadrons? The parton- 
parton interaction also diverges, and the infrared divergence in the two-body 
interaction exactly cancels the infrared divergence in the one-body operator 
for color singlet states. 

Of course, the cancellation of infrared divergences is not enough to obtain 
confinement. The cancellation is exact regardless of the relative motion of the 
partons in a color singlet state, and confinement requires a residual interac- 
tion. The 0(a) QCD hamiltonian contains a logarithmic potential in both 
longitudinal and transverse directions. There is no rigorous demonstration 
that the confining interaction is linear, and a logarithmic potential is of inter- 
est phenomenologically for heavy quark bound states.Ea I would be delighted 
if a better light-front calculation produces a linear potential, but this may not 
be necessary even for successful light hadron calculations. 

The calculation of how the quark self-energy changes when a similarity 
transformation lowers the cutoff on energy transfer is almost identical to the 
electron self-energy calculation. We find the one-body operator required by 
coupling coherence, 



where C F = {N 2 - l)/(2N) for a SU(N) gauge theory. 

The calculation of the quark-antiquark interaction required by coupling 
coherence is also nearly identical to the QED calculation. Just as in QED 
the coupling coherent interaction induced by gluon exchange above the cutoff 




(45) 
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partially cancels instantaneous gluon exchange. For the discussion of con- 
finement the part of V co h that remains is not important, because it produces 
the short range part of the Coulomb interaction. However, the part of the 
instantaneous interaction that is not canceled is 

/ >. 2 



(46) 



Vinstant = ~8g A C F \f Pi vt kjkj y— J 5 ai<J2 S Xl x 2 
x 9(\p+ - p+\ - eP+) 



x exp 



-2s(Pi -P 2 ~Q )(k 2 ~K -q 



Note that this interaction contains a cutoff that projects onto exchange 
energies below the cutoff, because the interaction has been screened by gluon 
exchange above the cutoffs. This interaction can become important at long 
distances, if parton exchange below the cutoff is dynamically suppressed. In 
QED I argued that this singular long range interaction is exactly canceled by 
photon exchange below the cutoff, because such exchange is not suppressed 
no matter how low the cutoff becomes. Photons are massless and experience 
no significant interactions, so they are exchanged to arbitrarily low energies 
as effectively free photons. This cannot be the case for gluons. 

For the discussion of confinement, place the most singular parts of the 
quark self-energy and the quark-antiquark interaction in Hq. To see that 
all infrared divergences cancel and that the residual long range interaction is 
logarithmic, study the matrix element of these operators for a quark-antiquark 
state, 

a\ rs 

<f>(p,a,r;k,\,s)b r Hp,a)d s i(k,\)\0) , (47) 

where r and s are color indices and <f> is a color singlet. 

When the expectation value of the hamiltonian is taken using this state 
there are divergences as e — > in both the expectation value of the self-energy 
and the surviving piece of instantaneous gluon exchange. These divergences 
cancel exactly for any color-singlet state. The cancellation resembles what 
happens in the Schwinger model. If the state is a color octet the divergences 
are both positive and cannot cancel. Since the cancellation occurs in the 
matrix element, we can let e — * before diagonalizing Hq. 

The fact that the divergences cancel exactly does not indicate that con- 
finement occurs. This requires the residual interactions to diverge at large 
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distances. It is easily shown that for large longitudinal separations the inter- 
action becomes: 



V(x-) 



g\c F K 2 



X 



I) 



(48) 



2tt 2 



At large transverse separations it becomes: 



V(x±) 



ln(|x. 



I) 



(49) 



The strength of the long-range logarithmic potential is not spherically sym- 
metrical in these coordinates, with the potential being larger in the transverse 
than in the longitudinal direction. Of course, there is no reason to demand 
that the potential be rotationally symmetric in these coordinates, because ro- 
tations are dynamical and are supposed to alter the Fock space composition 
of states in addition to rotating the state in a given Fock space sector. 

Had we computed the quark-gluon or gluon-gluon interaction, we would 
find essentially the same residual long range two-body interaction in every 
Fock space sector, although the strengths would differ because different color 
operators appear. In QCD gluons have a divergent self-energy and experience 
divergent long range interactions with other partons if we use coupling coher- 
ence. In this sense, the assumption that gluon exchange below some cutoff is 
suppressed is consistent with the hamiltonian that results from this assump- 
tion. To show that gluon exchange is suppressed when A — > Aqcd, rather 
than some other scale (i.e., zero as in QED), a non-perturbative calculation of 
gluon exchange is required. This same confinement mechanism appears using 
the similarity transformation developed by WalhoutEa 

I provide, pQnly a brief summary of our heavy quark bound state 
calculations ,ESEil and refer the reader to the original articles for details. We 
follow the strategy that has been successfully applied to QED, with modifica- 
tions suggested by the fact that gluops experience a confining interaction. 

For heavy quark bound states Eil we can simplify the hamiltonian by 
making a nonrelativistic reduction and solving a Schrodinger equation. We 
must then choose values for A, a, and M. These should be chosen differently 
for bottomonium and charmonium. The cutoff for which the constituent ap- 
proximation works well depends on the constituent mass, as in QED where 
it is obviously different for positronium and muonium. In order to fit the 
ground state and first two excited states of charmonium, we use A = 2.5GeV, 
a = 0.53, M c — 1.6GeV. In order to fit these states in bottomonium we use 
A = 4.9GeV, a = 0.4, and Mb = 4.8GeV. Violations of rotational invari- 
ance from the remaining parts of the potential are only about 10%, and we 
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expect corrections from higher Fock state components to be at least of this 
magnitude for the couplings we use. 

These calculations show that the approach is reasonable, but they are not 
yet very convincing. There are a host of additional calculations that must be 
done before the success of this approach can be judged. 
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